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Abstract 



Observational evidence for dark matter can be explained by Weakly Interacting Massive Particles (WIMPs). These dark matter 
particle candidates could indirectly be detected through the observation of signals produced as part of WIMP annihilations or 
decays. Latest results from indirect searches for WIMPs are reviewed. Current and planned experiments are presented and their 
prospects and discovery potential discussed. 
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1. Introduction 



Despite overwhelming evidence that it composes the 
vast majority of the mass in the Universe, dark matter's 
particle properties literally remain in the dark. Iden- 
tifying the mysterious nature of dark matter is one of 
today's most pressing scientific problems and is being 
sought for using colliders, direct-detection experiments, 
and powerful indirect detection techniques. WIMPs - 
Weakly Interacting Massive Particles (denoted x ) are 
attractive candidates for dark matter and naturally arise 
in many theories beyond the standard model of parti- 
cle physics, which were developed to explain the origin 
of electroweak symmetry breaking and solve the gauge 
hierarchy problem (TJ. 

This review is structured in the following way. Evi- 
dence for dark matter and expected signals are discussed 
first, followed by a review of experiments and recent re- 
sults. Prospects for the detection of WIMPs with ex- 
isting and proposed experiments are discussed before 
concluding. 

2. Evidence for dark matter and signals 

2. 1 . Evidence for dark matter 

Evidence for the existence of dark matter can be ob- 
tained at all scales from the motion of stars in dwarf 



spheroidal galaxies, galactic rotation curves, virial ve 
locities in clusters of galaxies, gravitationally lensed 
galaxies, baryon acoustic oscillations (BAO), to im- 
prints on the cosmic microwave background (CMB). 
N-body simulations of dark matter can reproduce ob- 
served large scale structures in our universe and can also 
be used to obtain expected dark matter distributions in 
Milky Way like galaxies. Dark matter only simulations 
yield spherically symmetric halo density profiles p(r), 
which describe the average dark matter density as func- 
tion of the distance r from the Galactic center. The 
impact of baryons, the behavior for small r (cusp-core 
problem), and the impact of sub- structure are still top- 
ics of debate. Figure[T]summarizes the observational ev- 
idence and lists astronomical objects that are expected 
to contain high dark matter densities and therefore are 
natural targets for indirect searches. 



2.2. Thermal relic 

If dark matter is a WIMP that is a thermal relic of 
the early Universe, then its total self-annihilation cross 
section averaged over the velocity distribution (cr^v) is 
revealed by the observed dark matter abundance. Evolu- 
tion of the number density n is determined by the com- 
petition between production and annihilation with stan- 
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Figure 1 . Observational evidence for dark matter and potential target 
that are expected contain significant amounts of dark matter. 



squared WIMP density over the integration solid angle 
AQ, defined as 

J(A£l)= \ dQ! \ p 2 (r(l,(f>'))dl(r,(f>') (3) 

JAQ(0,0) Jlos 

Note, that the dark matter annihilation cross sections 
could be velocity-dependent l3l . 

Expected fluxes of other cosmic messenger particles 
(y-rays, p/p, e ± , ...) can be obtained similarly, but often 
require extensive modeling of diffusion and propagation 
effects. WIMPs can be searched for with all these mes- 
sengers, but y-rays and neutrinos provide best detection 
prospects as spectral and directional information can be 
directly linked to the WIMP annihilations. Searches in 
the CMB are not discussed here ID. 

3. Instruments 



dard model particles. 



dn 

— +3Hn 
dt 



d(na 3 ) 
a 3 dt 



(cr A v) [n 2 eq - n 2 ) 



(1) 



The Hubble parameter, //, provides a measure of the 
universal expansion rate and n eq is the equilibrium abun- 
dance n eq = g x (mT/(27i)) 3/2 exp(-m/T). In equilibrium 
n decreases exponentially as the Universe expands and 
cools. When the^ abundance becomes very small, equi- 
librium can no longer be maintained and WIMPs freezes 
out [2], setting the natural scale for (cr A v>. 

2.3. Indirect Detection 

Indirect signals from WIMP annihilations depend on 
the square of the dark matter density, p, the annihilation 
products, and propagation effects. As an example we 
discuss the differential neutrino flux from WIMP anni- 
hilations of mass m x , which is given by: 



d<3> 

-=(£,0,0) 
dE 



1 <cr A v) 
4n 2m^ 



(2) 



where ^ is the differential neutrino multiplicity per an- 
nihilation weighted by the branching fraction Bf. This 
annihilation part encodes the particle physics in form 
of the energy spectrum. Intrinsic WIMP properties 
could be derived from spectral shapes, as well as the 
WIMP mass identified via end-point analyses. The "As- 
trophysics factor" / has a characteristic spacial depen- 
dence. / is the integral over the line of sight (los) of the 



PAMELA - Payload for Anti-Matter Exploration 
and Light-nuclei Astrophysics, is a satellite-borne mag- 
netic spectrometer attached to the Resurs-DKl satel- 
lite. Charged particles passing through the instrument 
aperture are deflected by the field of a neodymion-iron- 
boron permanent magnet. The resulting curvature of 
charged particle trajectories is measured precisely by 
a tracking system consisting of 6 plane double-sided 
silicon micro-strip tracker and depends on the particle 
charge and magnetic rigidity (R = pc/\Z\e). The en- 
ergy and interaction topology of particles is measured 
with a sampling imaging calorimeter, in which pairs 
of orthogonal ministrip silicon sensor planes are inter- 
leaved with tungsten absorber plates. The main purpose 
of the calorimeter is to distinguish e ± from p/p and He. 
Plastic scintillation counters act as an anti-coincidence 
system to reject particles that have entered through the 
sides of the instrument. Further detector components 
are a time of flight system, neutron calorimeter, and bot- 
tom scintillator l5l . 

Fermi-LAT - The Fermi Large Area Telescope is a 
pair-conversion telescope and the primary instrument 
on the Fermi Gamma-ray Space Telescope launched on 
June 11, 2008. It consists of three detector subsystems: 
(1) A tracker/converter with 18 layers of paired sili- 
con strips detectors interleaved with tungsten foils, (2) 
An 8 layer CsI(Ti) scintillation crystals calorimeter, (3) 
plastic scintillator tiles and wavelength- shifting fibers as 
anti-coincidence detector. Fermi-LAT is sensitive in the 
energy range of 20 MeV - 300 GeV. Fermi-LAT obser- 
vations started in August 2008 and first analyses were 
performed with the Pass 6 data event analysis scheme 
designed prior to launch. Pass 7 has been improved 
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by accounting for known on-orbit effects. Pass 8 is 
still being developed, but it is redoing everything from 
the event reconstruction up l6l . 

IACTs - Imaging Atmospheric Cherenkov Tele- 
scopes (IACTs) Q detect the Cherenkov light emit- 
ted in atmospheric particle showers with duty cycles of 
about 15% at a field of view (FOV) between 3.5° to 5.0° 
depending on the instrument. By using the atmosphere 
as a target their effective area is about a factor 500 larger 
than that of Fermi-LAT (~ lm 2 ). VERITAS, opera- 
tional since September 2007, is located at the Whipple 
Observatory Site (Southern Arizona) at an altitude of 
1250 m. It provides sensitivity to a wide range of ener- 
gies (150 GeV - 30 TeV) through stereoscopic imaging 
with four 12-meter telescopes. MAGIC |8], located at 
the Canary Island of La Palma at 2200 m, provides a low 
energy threshold of about 50 GeV with Stereo IACTs 
with two 17 m telescopes. Regular stereo observations 
were performed since Fall 2009. The High Energy 
Stereoscopic System (H.E.S.S.) is located at Khomas 
Highland of Namibia at an altitude of 1800 m consists 
of an array of four 13 -meter IACTs and achieves an en- 
ergy threshold of about 200 GeV. Recently, the 600m 2 
(28-meter-sized mirror) telescope H.E.S.S. II has been 
added. 

IceCube - The IceCube neutrino telescope instru- 
ments a volume of about one gigaton of Antarctic ice 
beneath the surface of the Geographic South Pole with 
5160 digital optical modules (DOMs) on 86 strings l9l . 
Each DOM contains a 10-inch PMT. IceCube exploits 
the good optical properties of the ice beneath the South 
Pole to detect neutrinos through the Cherenkov light 
emission of secondary particles produced by neutrino 
interactions. 

Construction was completed after 6 years in Decem- 
ber 2010. Data acquired during the construction phase 
has been analyzed. 

ANTARES - The ANTARES neutrino telescope (TO) 
is located between the depth of 2000 m - 2475 m in 
the Mediterranean Sea, roughly 40 km offshore from 
Toulon in France. The full detector consists of twelve 
vertical lines equipped with a total of 885 10-inch PMTs 
(R7081-20 from Hamamatsu). 

Super-K - The Super-Kamiokande detector ifTTIl is 
the world's largest underground water Cherenkov de- 
tector, located in Kamioka mine (1000 m of rock over- 
burden) in Hida-city, Gifu, Japan. The cylindrical 50 kt 
with a diameter of 39.3 m and a height of 41.4 m began 
operating in April 1996. Cherenkov light is used to re- 
construct the position of neutrino interactions in the de- 
tector (contained events) and up-going muons passing 
through the detector, which originate from muon neutri- 
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Figure 2. Positron fraction as function of the electron energy. 
PAMELA, HEAT, AMS-01, and Fermi-LAT data show a rising frac- 
tion above 10 GeV. The line indicates the expected positron fraction 
as obtained from standard electron-positron models. 

nos interacting outside the detector volume. 

4. Summary of current results 

4.1. Charged cosmic-rays 

WIMP annihilations in the Milky Way halo could 
yield positron, antiproton and antideuteron signals. The 
complex interplay of these messengers is discussed 
elsewhere (e.g. 1 12]) and we focus here on positrons. 

Measurements of the positron fraction cp(e + )/((p(e + ) + 
(p(e~)) show an unexplained rise above 10 GeV, that is 
inconsistent with standard secondary production mod- 
els Ifl3l . The excess, first observed by PAMELA fl4l 
following initial hints by HEAT and AMS-01, has been 
confirmed by Fermi-LAT by back tracking leptons in the 
Earth's magnetic field. The observed positron anomaly 
is shown in Fig. [2] (HI . 

Due to the small gyroradius of electrons in the 
Galactic magnetic fields, electrons lose energy quickly. 
Hence, the source of the lepton excess must be of lo- 
cal nature and originate within about one kpc of the 
Sun. It could be due to unaccounted or limitations in 
the modeling of nearby astrophysical sources (e.g. pul- 
sars fl6l . supernova remnants El), or be a first hint of a 
dark matter signal. If the positron excess in combination 
with the energy spectra of electrons and positrons above 
300 GeV measured by Fermi [El, H.E.S.S. [19], and 
ATIC l20l are interpreted as originating from WIMP an- 
nihilations or decays, then it indicates leptophilic mod- 
els that require enhancement in annihilation rate com- 
pared to the thermal relic cross section. Such scenarios 
are presently disfavored by y-ray, radio, and neutrino 
observations. 
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4. 2. Gamma-rays / radio 

Combined Fermi-LAT observations of various dwarf 
spheroidal galaxies have constrained WIMP self- 
annihilations with the thermal relic cross section for 
WIMP masses below 10 GeV I2T1L limits significantly 
weaken above on TeV. MAGIC has investigated dwarfs 
Willman 1 l22l and Segue 1 and obtained limits around 
10" 22 cm 3 s -1 and 10" 23 cm 3 s -1 for annihilation in yu + yu" 
and r + r", respectively, for WIMP masses between 
200 GeV and 1.2 TeV |53). VERITAS observed no sig- 
nificant gamma-ray excess from Segue 1 [24] and four 
dwarf spheroidals and derived an upper limit on the 
gamma-ray flux to constrain the self-annihilation cross 
section. The most stringent limit of 10" 22 cm 3 s _1 is ob- 



tained for m v = 300 GeV (25 



10" 21 cm 3 s 



H.E.S.S. reported cross 
1 to 10" 22 cm 3 s- 1 



section constraints from 
from the Sculptor and Carina dwarf galaxies depend- 
ing on the assumed halo model |26l . Limits can also 
be obtained from observations of y-rays, produced via 
final state radiation, inverse Compton scattering, or syn- 
chrotron radiation, of the Galactic Center and Galactic 
Ridge regions, as well as radio observation of the Galac- 
tic Center (27). Cosmic X-ray data further provides con- 
straints on WIMP annihilation l28l . 

Spectral features such as a spectral line can be a con- 
vincing signature for WIMP annihilations and could be 
produced in XX ~~ * 77 '» hy, and Zy. In a search by 
the Fermi-LAT collaboration in 2 yrs of Pass 6 data, 
no signal was observed in the ROI (region of inter- 
est), which consisted of the entire sky, but excluded the 
Galactic plane and known y-ray sources (1FGL). As- 
suming the background is described by a power law 
with a spectral index free to vary, limits were derived on 
(cr A v) l29l . The obtained limits are in mild tension with 
claims of an observation of a 130 GeV y-ray line near 
the Galactic center in 4 years of Pass 7 data l30ll3TTl . 
While it is unlikely that the line originates from astro- 
physical background l32ll , it could be an instrumental 
effect, for example as the result of a non-linear energy 
mapping. Fortunately, the energy mapping can be tested 
in Earth limb data. While the Fermi-LAT collaboration 
will have to clarify this topic, initial studies show indi- 
cations of non-linear effects l33l . H.E.S.S. II data or 
radio data could eventually also be used confirm or re- 
fute the line in independent measurements. For anni- 
hilations with final states Zy or hy, relativistic e ± are 
expected, which would generate synchrotron radiation 
when interacting with Galactic magnetic fields. Exist- 
ing radio data in the Galactic Center is already in minor 
tension with the presence for a contracted NFW profile. 
Currently running and future radio telescopes like Long 



Wavelength Array, LOFAR and SKA have great poten- 
tial in resolving this question further l34l[35l . 

4.3. Neutrinos 

The indirect search for WIMP-induced neutrinos 
aims to detect Galactic signals similar to the searches 
with gamma-rays, but also from self-annihilating 
WIMPs captured by the Sun and Earth. IceCube has 
searched for signals from the Galactic halo [36] and 
Galactic center l37l and improved upon theoretical pre- 
dictions l38l . Tight constraints were also derived from 
dwarf spheroidal galaxies | 39 ] and the Virgo cluster. 
Figure|3]shows a comparison of these present bounds on 
the dark matter self-annihilation cross section as func- 
tion of the WIMP mass for neutrinos and gamma-rays. 
Neutrinos are in particular competitive for m x > 1 TeV, 
with best sensitivity achieved by clusters of galaxies if 
substructure is taken into account l40ll . 
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Figure 3. Comparison of gamma-ray bounds with those from neu- 
trino telescopes and a dark matter model motivated by the PAMELA 
positron excess forxx ~~ * T+T ~ El- 

WIMPs could accumulate in the Sun or Earth and 
give rise to detectable neutrino signals. Energy loss in- 
duced by an initial scatter of a WIMP on a nucleon in the 
Sun can lead to the gravitational capture in the Sun. The 
probability for such an interaction, which is the same 
underlying physics process as being searched for in di- 
rect detection experiments, depends on the WIMP nu- 
cleon scattering cross section. WIMPs accumulate in 
the Sun and start annihilating. The annihilation rate de- 
pends on the amount of dark matter in the Sun. The 
annihilation rate steadily increases with the number of 
thermalized WIMPs near the center of the Sun up to a 
point where it becomes equal to the capture rate. At this 
point equilibrium has been reached and the annihilation 
rate is independent of the self-annihilation cross section. 
The neutrino flux from the Sun hence depends only on 
the capture rate, which can then be linked to the WIMP 
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nucleon scattering cross section. As the Sun is primary 
a proton target, in particular tight constraints can be 
derived on the spin-dependent WIMP-proton scattering 
cross section cr xp . Figure [4] shows IceCube's sensitiv- 
ity with one year of data collected with the 79-string 
detector. The dataset has been divided in three inde- 
pendent categories (summer, winter low-energy, winter 
high-energy) and is later combined. 
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Figure 4. IceCube 79-string detector solar WIMP sensitivity PTl . 

Super-K searched the up-going muon sample for sig- 
nals in 3109 days and derived limits l42l . that were now 
improved upon after including fully and partially con- 
tained events (see Fig. [5}, which are of importance for 
low-mass WIMPs 1431 . The preliminary limits from 
Super-K are compared to those from ANTARES and 
IceCube in Fig. [6] 

Limits on the WIMP-nucleon scattering cross section 
can also be deduced from limits on mono-jet and mono- 
photon signals at hadron colliders, however, they de- 
pend strongly on the choice of the underlying effective 
theory and mediator masses l44l . 
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Figure 6. Constraints on the spin-dependent WIMP-proton scattering 
cross section. 
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Figure 5. Super-K solar WIMP analysis: Distribution among fully 
contained (FC), partially contained (PC), and up-going muons (upmu) 
events expected at final analysis level for a WIMP of given mass an- 
nihilating into bb. For WIMPs masses below 50 GeV the dominant 
signal is expected from contained events (FC+PC) [45]. 



Figure 7. Sensitivity of a next generation neutrino detector assuming 
5 Megaton years of data | 43 1. 

A conclusive test of many low-mass dark matter sce- 
narios, a more precise study of atmospheric oscillation 
parameters, and an enhanced sensitivity towards super- 
nova burst neutrinos would require a very large neutrino 
detector with a low energy threshold. Such a detector 
could be constructed in two phases at the geographic 
South Pole, making use of the existing infrastructure, 
good optical properties of the naturally occurring detec- 
tor medium and support structure, and benefit from the 
IceCube detector to veto atmospheric muons. 

A mega-ton-sized ring-imaging detector could be 
constructed in two stages. The first stage (PINGU 
- Precision IceCube Next Generation Upgrade) would 
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consist of an upgrade to the IceCube-DeepCore detec- 
tor l46l using existing technology complemented with 
calibration devices and novel optical modules. Physics 
results are guaranteed by relying on proven IceCube 
sensors, while the performance of new technologies to- 
wards stage 2 can be evaluated. For stage 1 it is en- 
visioned to deploy about 20 string during two seasons. 
The aim is to achieve an energy threshold of a few GeV 
for this multi-mega-ton detector (47l[48]|. A stage 2 de- 
tector (MICA - Mega-ton Ice Cherenkov Array), con- 
sisting of the order of 100 strings, using a technology 
choice based on the performance of the stage 1 array 
would then aim at constructing a large ring-imaging de- 
tector. It is envisioned to reconstruct individual events 
above a threshold on the order of 100 MeV and use 
collective event information to detect supernova burst 
neutrinos. Primary physics motivations are dark mat- 
ter searches, neutrino oscillation studies, and increased 
sensitivity towards supernova burst neutrinos. Exten- 
sions reaching proton decay could possibly be contem- 
plated. Other physics topics include, but are not lim- 
ited to: Sterile neutrinos and neutrinos from Galactic 
sources. 

The deep ice at the geographic South Pole possesses 
good optical properties below 2100 m and a high radio- 
purity. The absorption length at 400 nm is about « 
155 m and effective scattering length is on the order of 
^scat ~ 47 m. Uranium ( 238 U) and Thorium ( 232 Th) con- 
taminations are very low at 10" 4 ppb and Potassium at 
0.1 ppb in the Antarctic ice l49l . The combination of 
low installation costs and the ability to build a contigu- 
ous detector, makes the South Pole an ideal site. How- 
ever, the maximum density of instrumentation is deter- 
mined by the installation procedure and will ultimately 
determine what photo coverage can be achieved. While 
the stage 1 detector can rely on the existing IceCube hot- 
water drilling technology for the stage 2 modifications 
are likely necessary. Drilling and deployment costs are 
expected to be below 10% of the total costs of the ar- 
ray, making the "excavation cost" component a moder- 
ate one. 

IceCube digital optical modules (DOMs) l50l are 
functioning extremely well, which is underlined by the 
fact that the number of DOMs that fail commissioning 
is at a percent level and the number of lost DOMs after 
successful freeze-in and commissioning is a fraction of 
a percent. The IceCube detector is operating very stable 
and shows detector uptimes of about 99%. DeepCore 
utilizes 252 mm diameter Hamamatsu R7081MOD (su- 
per bialkali photocathodes), which are identical to the 
standard IceCube PMTs (R7081-02), but with a quan- 
tum efficiency that is increased by 40% at X - 390 nm. 



While, the physics goals of the stage 1 detector are 
achievable with the existing DeepCore sensors, it is in- 
tended to utilize also new photon detection technolo- 
gies, with the goal to demonstrate the potential for re- 
constructing Cherenkov ring fragments. Developed for 
KM3NeT (SB, multi-PMT optical modules, could be 
adapted for the use in the ice. Other optical devices 
utilizing wavelength shifter techniques to increase the 
photo sensitive area in a cost effective manner are also 
under consideration. 



6. Future prospects and outlook 

6.1. Gamma-rays 

Fermi-LAT has produced already extremely compet- 
itive bounds on (cr A v), future analyses will benefit from 
larger statistics and potential discoveries of new ultra- 
faint dwarfs. Figure[8]shows a comparison of the Fermi- 
LAT dwarf spheroidal analysis lf2TTl compared to the 
expected sensitivity of the Cherenkov Telescope Array 
(CTA) (52), which is currently in the design phase. In 
the future GAMMA-400 (53) with an effective area of 
4m 2 , an angular resolution of ~ 0.01° at E 7 =100 GeV 
and an energy resolution of 1% could surpass Fermi- 
LAT. GAMMA-400 could be launched as early as 2018. 
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Figure 8. Comparison CTA sensitivity |52 | to current (2 years) Fermi- 
LAT exclusion limit |21| and scaled to the 10 years expectation. 

6.2. Neutrinos 

IceCube bounds presented at this conference only 
used the partially instrumented detector. More than 
a year of high-quality physics data has already been 
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Figure 9. Extra-galactic gamma-ray flux as observed by Fermi-LAT 
and predictions for neutrino fluxes testable by IceCube (54). 

collected with the full IceCube detector, including 7 
more strings (2 of them for DeepCore), further lower- 
ing its threshold. This data could also be used to test 
un-explained features in the extra-galactic gamma-ray 
fluxes observed by Fermi-LAT, that could hint at high 
mass dark matter l54l (see Fig. [9]). 

Exciting prospect for dark matter searches also ex- 
ist at next-generation neutrino detectors such as Hyper- 
Kamiokande [55 ], LENA (56), and PINGU. WIMP sce- 
narios motivated by DAMAs annual modulation sig- 
nal l57l and isospin- violating scenarios l58l motivated 
by DAMA and CoGeNT signals would be testable at (~ 
lGeV) - threshold neutrino detectors l43l (see Fig. [7]). 
Further new detection channels such as low-energy neu- 
trinos originating from hadronic particle showers cre- 
ated in WIMP annihilations in the Sun could enhances 
sensitivities and provide new ways to search for dark 
matter. 

6.3. Low-energy neutrinos 

The large existing Super-K dataset and that of pro- 
posed future low-energy threshold neutrino detectors is 
sensitive to a new WIMP detection channel | 59, 60). A 
large number of low-energy neutrinos are expected from 
any annihilation with hadronic components in its final 
states in the Sun. Hadrons interact in the dense solar 
medium, producing a hadronic shower, which will re- 
sult in a large number of pions (the pion multiplication 



effect is shown in Fig. 10 ). Positive pions and muons de- 
cay at rest, producing low-energy neutrinos with known 
spectra, including v e through neutrino mixing. A sig- 
nal could be detected by Super-K (see Fig. [TT]) and in 
the future provide a new probe of WIMP-nucleon scat- 
tering. Compared to other methods, the sensitivity is 
competitive and the uncertainties are complementary. 
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Figure 10. Low-energy neutrino yield (summed over flavors) per 
WIMP annihilation ( <\[s = 2m x ) in the solar core, obtained by sim- 
ulating pion-induced hadronic showers in the solar medium is com- 
pared to the high-energy neutrino yield for t + t~ (59). 
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Figure 1 1 . Example of a low-energy v e signal from solar WIMP an- 
nihilation in Super-K, along with the measured background (4.1 live- 
time years). The signal shape is independent of WIMP properties, 
and its normalization scales with o*j® (here chosen to be at the edge 
of exclusion) (59l . 

7. Conclusions 

Indirect searches have resulted in tight constraints 
on fundamental properties of WIMPs. Neutrinos are 
highly competitive for extended (Galactic halo/center, 
Clusters of Galaxies, ...) and point-like sources (Dwarf 
spheroidal, ...) for energies above about 1 TeV. Gamma- 
rays are more sensitive in the detection of lower-mass 
WIMPs. Any detection of WIMPs would likely require 
an independent observation with both messengers. 

WIMP annihilations in the Sun provide a discovery 
channel for neutrinos through a striking signature. Due 
to the opacity of the Sun to high-energy neutrinos orig- 
inating in the center, Solar WIMP signals are detected 
in the energy range below 100 GeV. Neutrino detectors 
(IceCube and Super-K) provide the world's best limits 
and their sensitivity continues to improve as more data 
is collected and new detection channels are investigated. 
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